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^ ■ Abstract 

(N ■ 

In June 2012 on a conference in Bielefeld, after the author made the presenta- 
tion of his theory of nonlinear Markov processes, Tom Kurtz asked him whether his 
methods would allow to get well-posedness for nonlinear McKean-Vlasov type dif- 
fusions with coefficients depending on the distribution via the median. The present 
paper is the answer to this question. The aim here is not to find the weakest possible 
conditions for the affirmative answer, but a sensible general result with a quick and 
transparent proof. 
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1 Main results 

By V we denote the gradient with respect to the spatial variable x G R d , by \\u\\l the 
t— i ■ norm of functions in Li(R d ), by S x the Dirac measure at x G R d , by \v\ the Euclidean 
j> ! norm of a vector v G R d . 

For a positive integrable function u on R and an a G (0, \\u\\l) the quantile Q a (u) is 

$_i ; denned as the number satisfying 

& '. 

u(x) dx = a. 



- 1—1 



If a = 1/2, the value Q a (u) is called the median of u. In finances, the values of Q a (u) for 
small a specify the so-called Value at Risk (VaR). Other particular cases include quartiles 
and sextiles. 

More generally, for any positive function u on R d with = 1 and any vector a = 
(ai, • • • , a.d) with all acj G (0, 1), let us define the quantile vector Q a (u) = (Q\, • • • , Qt)( u ) 
with the coordinates uniquely specified by the equations 

u(x) dx = atj, j — 1, • • • , d. (1) 

{x&K d :x j <Q j a {u)} 
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The problem addressed in this paper concerns the well-posedness of the nonlinear SDE 

dX t = b(t, X t , Q a (C(X t ))) dt + <r(t, X t , Q a (C(X t ))) dW t , (2) 

where Wt is a standard d- dimensional Brownian motion, b, a are, respectively, vector- 
valued and (dxd square symmetric)-matrix- valued continuous functions on R + x R d x H d , 
and C(X) denotes the distribution density of a random variable X. Such SDE arise in 
the analysis of interacting particle systems related to asset pricing evaluation, see [3J. 
As an auxiliary problem we shall use the equation 

dX t = bit, X t , u(t)) dt + a{t, X t , u(t)) dW t , (3) 

with a given continuous curve u(t) in R d . By Ito's formula, the corresponding diffusion 
has the generator 

L t = ^a 2 (t, x, + b{t, x, (4) 



in one-dimensional case and 



1 d 2 d 

L t ^ - ^aijUk^x^it))^-^ + ^b j {t,x,u{t))—, (5) 



or, in vector notations, 

L t = i(a 2 (t, x, w(t))V, V) + (b(t, x, u(t)), V) (6) 

in case of arbitrary d. 

Recall that a square matrix-valued function a is called uniformly elliptic if 

m^iei^MM^o^Hei 2 (7) 

for all (,i6 R d , t > and some constant m > 0. 

It is well known (and follows from the existence of a regular Green function, see below) 
that, for uniformly elliptic coefficients, the distributions of solutions to equation (j3J have 
densities u t with respect to Lebesgue measure (even if the initial distribution does not 
have a density), and the evolution of densities is known to satisfy the dual equation 

^-utix) = L* t u(x) = X - V V(t, x, u{t))u t {x)\ - V(6(t, x, u(t)),u t (x)). (8) 
The equation for the distribution densities of processes solving (j2J) on [0, T] becomes 

^u t (x) = ^V 2 [a 2 it, x, Q a (u t ))u t (x)\ - V(6(t, x, Q a (u t )),u t (x)), t e (0, T\. (9) 

The well-posedness for this functional PDE is equivalent to the well-posedness of 
nonlinear SDE ([2]) (In probabilistic language, solutions to (EJ) supply weak solutions to 
SDE but once the unique Q a (ut) is plugged into ([2]), its strong, path- wise, well- 
posedness becomes straightforward). 

It will be convenient to distinguish two cases, with a depending and not depending 
on u, the second case requiring stronger regularity assumptions on the initial values Uq. 
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Theorem 1.1. Suppose a 2 is uniformly elliptic, twice continuously differentiable in x 
(with bounded derivatives) and does not depend explicitly on u, b is a bounded function 
continuously differentiable in x (with bounded derivative), a,b are continuous functions 
of t and b is Lipschitz continuous function of u, with the Lipschitz constant k. Then, 
for any T > and any (strictly) positive u G C(R d ) fl Lx(R d ) there exists a unique 
classical solution u t G C(R d ) fl Li(R d ) of equation (JH]) with initial condition Uq, that is, 
it is continuous in t G [0, T] in the topology of Li(R d ), satisfies ([9]) fort G (0, T] and 
coincides with u for t = 0. 

For the second result we need Sobolev spaces H\ (R d ) defined as the spaces of integrable 
functions on R d with all derivatives up to and including order k being well defined in 
distribution sense and being again integrable functions with the norm ||«||#fc(Rd) defined 
as the sum of Lx-norms of u and all its partial derivatives up to and including order k. 

Theorem 1.2. Assume the same condition as above, but with a being also a Lipschitz 
continuous function of u (with the Lipschitz constant x). Assume also that a, b are twice 
continuously differentiable in x (with bounded derivatives). Then, for any T > and 
any (strictly) positive Uq G Hf(R d ) fl C(R d ), there exists a unique classical solution u t G 
C(R, d ) n Hf(R d ) of equation ([9]) with initial condition uq. 

Theorem 1.3. Under the conditions of either of Theorems li.il or 1 1.21 for any £ G H d , 

there exists a unique global (i.e. defined for all positive times) strong solution of SDE (T3]) 
with the initial condition X = £. 

Our method is robust in the sense that it can be applied to other processes, not 
necessarily diffusions. Let us illustrate this claim by applying it to stable-like processes. 
For simplicity we shall discuss the processes with a fixed stability index a G (1,2) and 
only with drift depending on a quantile. 

Let us work directly with PDEs and generators. The corresponding stable SDEs can 
be written using the technique from [2]. Let us consider the processes in R d governed by 
the time-dependent generator 

L t = (b(t,x,u(t)),V) - a(x)|A| Q / 2 . (10) 

The evolution of densities of the corresponding process are generated by the dual operator, 
that is they solve the equation 

^u t (x) = L*u(x) = \A\ a / 2 [a(x)u(x)} - V[6(t, x, u(t))u t (x)\. (11) 
Our goal is the evolution depending on the quantile, that is, the equation 

^u t {x) = \A\ a / 2 [a(x)u(x)} - V[b(t,x,Q a (u t ))ut(x)]. (12) 

Theorem 1.4. Assume that a(x) is a positive function, bounded from below and above 
and twice continuously differentiable (with bounded derivatives), and the vector-valued 
function b is twice continuously differentiable in x (with bounded derivatives), continuous 
in t and Lipschitz continuous in u>. Then, for any T > 0, any (strictly) positive Uq G 
Lx(R d ) n C(R d ) and any £ G H d , there exist unique solutions u t G C(R d ) fl Li(R d ) of 
equation with initial conditions Uq or 5^ . 



3 



2 Discussion 



The regularity of initial conditions and coefficients in Theorem 11.21 can definitely be weak- 
ened by using more advanced regularizing properties of heat equations. For instance, using 
the smoothing property for diffusion propagators (resolving operators) in the form 

l|tf*'7llfl?(R<) < C{t - S )- 1/2 ||/||^ ( R,) (13) 

that should hold under certain regularity assumptions for coefficients (though the author 
did not find a proper reference), one can directly improve Theorem 11.21 by allowing w to 
belong to Hl(R d ) (and not only Hf(R d )). 

Theorem 11.31 can be extended straightforwardly to a weakly dense (in the set of prob- 
ability measures) class of initial conditions having the property that their quantile is a 
single point. Say, if d = 1, on can take the initial distribution in the form Yl^iPj^iji as 
long as there exists k such that 

fc-l k 

S/'A < « < ^pA- 

3=1 3=1 

Then the quantile Q a (ut) will tend to as t — > 0. An interesting question remains on 
what can be said about processes with initial distributions with the median not uniquely 
defined, as say (5 a + 5b)/ 2 with any a ^ b. 

Theorem 11.41 can be extended in several directions. It has a straightforward extension 
to processes with variable stability index a(x) with corresponding required properties of 
transition probabilities following by (straightforward time-dependent extension of) Theo- 
rem 5.1 of [7] (see also Theorem 7.5.1 from [H]), as well as stable processes with variable 
spectral measure and/or perturbed by compound Poisson processes (say, tempered stable 
processes) with required properties of transition probabilities following by Theorem 4. 1 of 
[7J. By developing appropriate smoothing results for time dependent stable-like processes, 
various 'stable' analogs of Theorem 11.21 can be obtained. 

Our arguments are inspired by the approach developed systematically in [ID] , though 
technical details are problem-specific. The main impetus for this development was an at- 
tempt to get a general understanding of the appearance of nonlinear measure- valued evo- 
lutions as limits of Markov models of interacting particles (extending the well developed 
theory arising from the standard McKean-Vlasov diffusions). Recently, there appeared 
many publications combining mean-field interacting particle systems with control, see [T], 
[I], [S], [S], [12]; [32] for various directions. From the results of the present paper one can 
expect an interesting development of mean-field control theory with peculiar (but from 
practical point of view quite natural) dependence of control on empirical measures via 
such characteristics as median or VaR. 

The rest of the paper is devoted to the proof of Theorems 1.1-1.3. 
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3 Two-sided estimates for heat kernels and their con- 



sequences 

Let us denote by G a the heat kernel of the standard heat equation in R d with a diffusion 
coefficient a > 0, that is 

G a (t, x) = (2nta 2 )- d / 2 exp{-^}. (f 4) 

Recall that J G a (t, x) dx = 1 for all a, t > 0. 
For the general heat equation 

du I 

— = L t u, L t u(x) = -{a(t, x)V, V)u(x) + (b(t, x), Vu(x)) + c(t, x)u(x). (f5) 

the following fact is well known (see e.g. [15], [H] and references therein). 

Proposition 3.1. (i) Suppose a is uniformly elliptic, that is ([7]) holds, and continuously 
differentiable in x and a, b, c are bounded and continuous with respect to all their variables, 
so that 

sup max( | Va(t,x) |, \b(t,x)\, \c(t,x)\) < M, 

t,x 

then there exist constants o~i,Ci, i = 1,2, depending only on m,M,T such that the Green 
function (or heat kernel) of equation (fl~5|) (i.e. its solution G(t,x,s,l;) with the initial 
condition 5^ at time s) is well defined and satisfies the following two sided bounds 

C 1 G ai {t-s,x-£)<G(t,x,s,£,)<C 2 G a2 (t-s,x-£,), 0<s,t<T. (16) 

(ii) Moreover, if a is twice continuously differentiable in x and b,c are continuously 
differentiable (with al derivatives bounded), then G(t,x, s,£) is differentiable in x and 



d 

—G(t,x,s,0 



< Cr 1/2 G(t, x, s, 0, < s, t < T, (17) 



where the constant C depends only on the bounds for the derivatives and m, M, T. 

Remark 1. The known proofs of two-sided estimates (especially the lower one) are rather 
technical. Thus it is useful to note (especially for possible generalizations) that we need 
only a weaker local lower bound (which follows from small time and near diagonal asymp- 
totics of the Green function, see e.g. ' t 8]). Namely, for Proposition (i) we need the 



l.h.s. inequality fTl6|) for \x — £| < cy/t — s for at least one c > 0, and for Proposition \3.4 
(ii) we need it for \x — £| < cyjt — s for sufficiently large c > 0. 

Notice that the Green function for equation ([8]) with L* t is obtained by changing places 
of x, £ and t, s in the Green function of the corresponding equation with L t , which is of 
form (fl5l) . 

Let us deduce two corollaries, one from the upper bound and another from the lower 
bound of inequality (TTB"]) . Let us denote by U m: M,t(u ) the set of solutions u t (x), t > 0, 
of the Cauchy problems of all equations f|T5|) satisfying the conditions of Proposition 13.11 
with given m, M and a given initial condition m . Also set 

U m ,M,(0,T]{Uo) = Ut£(0,T]U mi M,t( U o) ■ 
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Proposition 3.2. For any uq G Li(R d ) and e > there exists K > such that 

\u(x)\dx<e (18) 

\x\>K 

for all u G U m! M,[o,T\(uo) (that is, the set of measures {u(x)dx}, u G U mj M,[o,T](uo) is 
tight). 

Proof. Let K be chosen in such a way that 

/ |^o(^)| dx < e, 

J\x\>K 

and let K be another positive constant. 

Recall that the solution to the Cauchy problem for equation (fT5l) can be written in 
terms of its Green function as 



u t (x) = J G(t,x,0,OMOd£. 
From the upper bound in f|T6|) we get 

/ u t (x) dx < C 2 / / G a2 (t, x — £)u (!;) d£dx 

J\x\>K+K J\x\>K+K Jn d 

<C 2 G a2 (t, y)u (£) d£dy + C 2 G a2 (t, y)u (£) d£dy, 

J\Z\>K,ye-R d J\y\>K,\(\<K 

which (by changing the variable of integration in the second integral to z = y/Vi) is 
estimated by 

C 2 e+ \\uq\\ l [ (27ra 2 2 )- d / 2 exp{-^}dz. 

J\z\>K/VT Z(J 2 

Clearly the r.h.s. can be made arbitrary small by choosing e > small enough and K 
large enough. □ 

Proposition 3.3. For any strictly positive continuous function uq on H d and any K > 

there exists 5 > such that 

u(x) > 5 (19) 
for all u G t/ m ,A/,[o,T](wo) whenever \x\ < K . In other words 

M{u{x) :\x\<K,uE U miM ,[o,T}(u )} > 0. 

Proof. From positivity and continuity of u we can find, for any K > 0, a number 5 > 
such that u (x) > 5 whenever \x\ < K + y/T. 

Next, from the lower bound in ([To]) (even its local version) we get 



u t (x) > Ci / G ai (t, x - u o{0 d£. 

J\x-£\<Vi 
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If |x| < K we have |£| < K + yT in the integral. Consequently, for \x\ < K, 

u t (x) > 8C1 / G ai (t,y) dy 
J\ y \<Vt 

> 5C 1 (2vrt ( x 2 )- d / 2 exp{--i I } / dy = dC^a 2 )^' 2 exp{~}V d , 

2cT i J\y\<Vt 2a i 

where Vd is the volume of the unit ball in R d . This lower bound proves the proposition. □ 
Combining these two facts we arrive at the following main result of this section. 

Proposition 3.4. Consider equation (|T5|) under the assumptions of Proposition 1X71 (i) 
and the additional condition that the solution preserves the L\-norm (like equation (jSJ) 
that we are mostly interested in). Let also arbitrary a = (ai, • • ■ , ay) with ay G (0, 1) and 
T > be given. 

(i) For any u G C(R d ) U Li(R d ) with \\uq\\l = 1, there exist K, 8, e > such that, for 
any u G U mM ,(o,T]{uo) , 

M{u(x) : max \xj\ < 2K} > 5, (20) 

3 

mm\Qi{u)\<K (21) 

3 

and 

I \u(x)\dx<e. (22) 

J{x£~R d : minj \xj\>K} 

The set of functions satisfying conditions (1201) - (1221 zs convex. 

(ii) There exist constants K, 5 > suc/i i/iai, /or any £ = (£i, • • • , £d) £ R d 
w G U mt M,t{8^), t G (0, T], 

inf{u(a;) : max \ Xj - &| < 2K} > 5r d/2 , (23) 

3 

min\QU8s) - t 3 \ < Kyft (24) 

3 

and 

i \u{x)\dx<e. (25) 

J {x&n, d :mmj \xj-$j\>Ky/i} 

Proof, (i) Recall that Q a are defined by ([I]). Let us pick up an e such that 

e < min(a!i, • • ■ , a d , 1 - ai, ■ ■ • , 1 - a d ). 

Applying Proposition E2] we can find K such that ( ITSl) or ( 1221) holds for all -u G £7m,M,[o,T]( u o) 
implying (12T]) . Applying Proposition 13.31 yields ( )20|) for some 5 > 0. The last statement 
is obvious. 

(ii) Notice that u G U mj M,t(^) is i n fact the heat kernel 

u = G(t, x, 0, £) 

for the diffusion equation of a considered class. Using the upper bound in (TIB]) , the free 
heat kernel ( 1141) and a well known (and easy to prove) estimate 

r°° r 2 K 2 

(2^r^J K exp{-_ }i<Cexp { -_ } 

for t G (0,T] with some constant C depending on T and a only, we obtain (125]) for large 
enough JC. This implies (I24p . The lower bound in ( 1T6|) implies (123]) . □ 
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4 Sensitivity analysis for diffusion equations 

Here we shall exploit the smoothness of the heat kernels of diffusion equations to prove 
Lipschitz continuity, in Li-norm, of the solutions to diffusion equations with respect to 
functional parameters. The results may be known, but the author did not find an appro- 
priate reference. 

Proposition 4.1. Consider two equations ffl5|) . specified by two families of operators 
L\,L\ with the coefficients ai,bi,cx and a 2 ,b 2 ,c 2 respectively, each satisfying the assump- 
tions of Proposition lff.il 

(i) Assume that the second order part of Lj,L 2 coincide, that is a% = a 2 . Then for 
any u$ G Li(R d ), the solutions u\,u 2 of the corresponding Cauchy problems satisfy the 
estimate 



u 



)-u 2 t \\ L < CVtsup(|&i(t,x) -b 2 (t,x)\ + \ Cl {t,x)-c 2 {t,x)\)\\u \\ L . (26) 



t.x 



(ii) Not assuming that a\ = a 2 , but instead requiring that all coefficients ai,bi,Ci be 
twice continuously differentiate in x (with bounded derivatives), it follows that for any 
uq G Hf(R, d ) the solutions ul,u 2 of the corresponding Cauchy problems satisfy the estimate 

\\u]-u 2 t \\ L < CVtsup(\ai(t,x)—a 2 (t, x)\+\bi(t, x)— b 2 (t, x)\+\ci(t, x)—c 2 (t, ^)|)||Mo||i?2( R d). 

t,x 

(27) 

Proof, (i) We apply the usual method for comparing semigroups or propagators. Denoting 
by U\' s the propagator solving equation flT5|) for L l , that is 



U!' s f{x) = J Gi(t, x, s, 0f(0d£, i = 1, 2, 



with Gi the heat kernel of the diffusion generated by Li, we see that they satisfy the 
identity 

Ul>° - V? = f ^-U l 2 s Uf ds = f U?{L\ - Ll)Uf. (28) 
Jo ds Jo 

From ffl7|) it follows that, for any u G Lx(R d ), 

||Vf/r'X||L<C S - 1/2 ||n || L 

with a constant C. Hence, from ff28l) and using the assumption that the difference L\ — L 2 S 
is a first order operator we conclude that 



\\{U^ -U^)u \\ L <C\\u4 L / s- 1 / 2 ds = 2Ct 1 / 2 \\u \\ L , 



as required. 

(ii) It follows by the same identity ( 1281) combined with the standard fact that the 
propagators (resolving operators) of uniformly elliptic diffusion equations are bounded 
operators in Hf(R d ) (here the assumption of the existence of the second derivatives of 
the coefficients is needed), and with the observation that 

\(Ll-I%)f(x)\ L < supila^t, x)-a 2 (t, x^ + lhit, x)-b 2 (t, x^ + lc^t, x)-c 2 (t, x)\)\\f\\ H 2 {Ild) . 

t,x 

□ 
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5 Proof of the Theorems 

Proof of Theorems \1.1\ and \1.2[ 

Let d = 1. Let us pick up a uq and then choose K,S,e > from Proposition 13.41 
Let ui,U2 G t/ mi Af,[o,T](wo)- Then the function u{K) = Q a (ui + h{u2 — Mi)) satisfies the 
equation 

[ui(x) + h{u2(x) — u\(x))] dx = a. 



Differentiating with respect to h yields 

u'(h) = —[ui(u(h)) + h(u2(u!(h))—ui(u(h)))}^ 1 / (u2{x) — u\{x)) dx . 



By Proposition 13.41 we deduce the following basic estimate 

f 1 1 

\Q a {u 2 ) - Q a {u x )\ = HI) - w(0)| = I J co'(h) dh\ < -\\ Ul - u 2 \\ L . (29) 
Similarly, for arbitrary d, we get for (Jj(h) = Q 3 a (ui + h(u 2 — Mi)) the equation 
uj'jih) / (wi + Zi^-wi))^)!^.^.^) JTcfo fe = / (w 2 (x) -mi(x)) dx, (30) 

implying, by ( l20|) and (EE]), the estimate 

l<£M - < ^=lll«i " Ml (31) 

for each j. 

Next, let C([0,T],R d ) denote the Banach space of -Revalued continuous functions 
on [0,T] with the usual norm \\u). \\ = sup t \u t \, and let C(K) be the cube of side 2K 
centered at the origin. For any u G C(K) let C wo ([0, T], C(K)) denote the convex subset 
of C([0, T], R d ) consisting of curves with oj given and with values in C(K). Let 

W = Q Q (U ). 

For a given curve co, G C Wo ([0, T], C(K)) let mJo;.](x) denote the solution at time t of 
equation ([8]) with the initial data uq(x). Let us define 

&t = Q a {u t [u]). 

By the strong continuity of the propagators solving equation ([H]) in Li and by (l3"Tj) we 
conclude that u)j depends continuously on t. Consequently, applying Proposition 13.41 we 
deduce that the mapping u, oj. is a mapping from C Wo ([0, T],C(K)) to itself. It is clear 
that the functions u t solve equation flU]) if and only if uj t = Q a (ut) is a fixed point of this 
mapping. Thus well posedness of equation ([HD is reduced to the problem of uniqueness 
and existence of this fixed point. 

Let uj x ,uj 2 be two curves in C Uo ([0,T],C(K)) and w 1 ,^ 2 their respective images. By 

(ED- 

\u\ - uj 2 t \ < jj^Wutlu 1 ] - u t \uj 2 }\\ L . 
9 



Moreover, from Proposition I4.1[ we get under the assumptions of Theorems 11.11 or 11.21 



the estimates 



\u t [u l ] -u t [u 2 }\\ L < VtCx sup \ul -u 2 \\\u \\l 

te[o,T] 



or 

\\utlu 1 ] -u t [u 2 ]\\ L < VtCx sup \u] - u 2 \\\u II u2 (Rd) 

te[o,T] 1 

respectively. In both cases this implies 



< v^C sup |w* -w s 2 | (32) 



se(o,i] 



with another constant C (depending on u ). Consequently, the mapping u, i->- a), is a 
contraction for small enough T. This implies the required well-posedness for small T. 
Global result follows by the usual iteration procedure completing the proof of Theorems 
O and PI 

Proof of Theorem \1.3[ Since the heat kernel (solution with the Dirac initial data) 
G(t, x,0,£) is a positive continuous function for all t > 0, and by Theorems 11.11 and 11.21 
it is enough to show the required well-posedness only for small times. Next, by ( 1301) . ( 13~TT) 
and the estimates (l2"3l. (EH), we find 



\Qi( u t) - ^K 1 )! < j^rTt d/2 t- {d - 1)/2 \\ul -u 2 \\l< J^TT^l - u 2 \\ L (33) 

for each coordinate j and 

uj(a:) = Gi(t,x, 0,0, « = 1, 2, 

the heat kernels of two equations jSJ) from the class specified by Proposition 13.11 (i). 

Next, under the assumptions of Theorem II .![ let us follow the same strategy by looking 
at the fixed points of the mapping w \—y w. from C^([0, T], C^{KyT)) to itself, where 
C^(Ky/T) is the cube of radius 2Ky/T centered at £. As in Proposition 14.11 (and noting 
that the norm of <5g is one) we get for heat kernels arising from a; 1 , a; 2 the estimate 

Hm^o; 1 ] — m 4 [w 2 ]||l < \ftC*c sup \u\ — Co> 2 |. 

se[o,t] 

Combining this with the previous estimate leads again to ( 132]) implying the contraction 
property of the mapping oj i— y Cj and hence the required well posedness. 

Let us turn to the assumptions of Theorem II .21 It will be more convenient now to 
work not with the quantiles u t , but with their transformations 

rh = t- 1/2 {0H-t)- 
Let us look for a fixed point of the mapping 

T). i — y to. i — y u. i — y cD, i — y fj. 

in the set C((0,T],C(K)) of continuous functions (0, T] !->■ C(K). To see that 77 1 — >- 77 is 
in fact a mapping in this space we note that fjt £ C(K) for all t (by Proposition 13.41 (ii)) 
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and that fjt is continuous for t > (by the proof of Theorem II . 3[) . It remains to obtain 
the contraction property. 

As in the proof of Proposition 14.11 (ii) (and using the fact that the second derivative 
of G(t,x, 0, £) with respect to x has singularity of order 1/t) we obtain 



\ut[uj l ] - uAuj 2 



r* i 

k 2 ]|U < Ck / -\uj\-ul\ds 
Jo s 



and hence 



\u t [co 1 (ri 1 )}-ut[io 2 (ri 2 }\\ L < CnVt sup \r, l s - rf s \. 

s£(0,t] 



Combining with f l3"3"j) we get 



\vl-V 2 \<ViC sup \vl-Vtl (34) 
se(o,i] 



which again implies the required contraction property. 



Proof of Theorem TJf. 



It goes along the same line. We just need the smoothing properties of equation ( Till . 
They are developed in [7j, [8] and j9], see also Chapter 7 of [UJ. It is shown (see Theorems 
3.1, 3.2 of [7j for time-homogeneous equations and Theorems 4.1, 4.2 of [9] for a straight- 
forward extension to time non-homogeneous case) that the transition probabilities (Green 
functions) G a of processes generated by operators ffTUj) satisfy the two sided estimate 

CiS a (t — s,x — £; a, b) < G a (t, x, s, £) < C^S^t — s, x — £; a, b), 0<s,t<T (35) 

with constants C\, C2 and S a denoting the stable density, i.e. the transition probability of 
the process generated by ffTUl) with constant a, b (Ci, C 2 ,a,b depend only on the bounds 
for a(x), a _1 (x), b(t, x, u) and their derivatives in x), and the derivative of G a with respect 
to x (or £) has the bound 

— -(t,x,s,0 < C(t- s)- 1/a GJt,x,s,0, 0<s,t<T (36) 
ax 

(see formula (58) of [9]). 

With these estimates and the standard estimates for S a (implying that S a is of order 
f-d/a f Qr a -fo^ Q £ rac [ ms ^i/a arounc l the diagonal, see e.g. books [7] or [UJ) the proof of 
Theorem 11.41 and all intermediate results repeats literally the proof of Theorem 11.11 
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